Entropy driven formation of Smectic C in system of Zig-zag shaped molecules 
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We have carried out Monte Carlo simulations of zig-zag shaped molecules where the the molecules 
is composed of three rigidly linked hard spherocylinder arranged in a zigzag fashion. By varying 
zigzag angle we have mapped out the whole phase diagram as a function of pressure, molecular 
density and zigzag angle fy. For ^ between 35° and 80° our model simulation exhibits SmC phase. 
This is the first conclusive evidence where steric interaction arising out of molecular shape alone 
induce the occurrence of SmC phase for a wide range of zigzag angle. For smaller *]/, a transition 
from tilted crystal (XT) to crystal is observed. 
PACS numbers: 61.30.-v, 61.30. Cz 



Liquid crystal phase is very sensitive to the 
molecular shape. Relating the macroscopic prop- 
erties of liquid crystal phase to the microscopic 
structures of the molecules is very complex prob- 
lem. On the other hand it is essential to under- 
stand this structure-properties relationship both 
from fundamental point of view as well as for 
different technological applications. Theoretically 
it is a very difficult task since it involves hun- 
dreds of molecules and their mutual interactions. 
Computer simulations has emerged as an impor- 
tant tool to investigate the dependence of the 
liquid crystal phase behavior to the structure of 
the constituent molecules. Various models has 
been used such as simple spherocylinder or hard- 
sphere chains which interact through hard or soft 
excluded- volume repulsion 0- 0. Q . molecules with 
ellipsoid shape interacting through Gay-Berne po- 
tential 0, , simple "bead-spring" representation 
of molecules interacting through Lcnnard-Jones 
potential • There are also simulations involving 
atomic level models of real liquid crystal molecules 

HQ. 

There are a number of theoretical works which 
gives a microscopic theories for the formation of 
smectic C phase [13 H E US El H3 • In McMil- 
lan theory it is the electric dipole-dipole interac- 
tions which produce tilt in the smectic C. Later 
quadrupolc-quadrupolc interaction was introduced 
as a perturbation of smectic A phase to get a smec- 
tic A to smectic C transition 20] . Motivated by the 
large discrepancy observed between X-ray and the 
optical data for the molecular tilt angle in smectic 
C materials, Bartolino eraJ. 1131 . resuming an idea 
from Guillon and Skoulios [lJI, proposed a molec- 
ular model (known as the zigzag model) derived 
from lyotropic systems: the rigid central core, op- 
tically anisotropic, imposes the tilt of the optical 
axis and the melted aliphatic end chains are on av- 
erage closer to the normal to the layers. Wulf [22^ 
while studying the zigzag shaped molecules con- 
cluded that the tilted smectic C is formed due to 
the steric interactions arising out of the packing 
arrangement of such molecules. There are num- 
ber of simulations concerning the occurrence of 
smectic C phase, but none of them are conclu- 



of model molecules comprising three rigidly linked 
Gay-Berne (GB) sites arranged in a zigzag fashion 
|T2l | and but did not find the occurrence of smectic 
C phase. Evidence of smectic C phase with vary- 
ing tilt angle from layer to layer was reported in 
simulation of zigzag shaped molecules comprising 
seven rigidly bonded soft spheres ^3] ■ One of the 
major goal of this paper is to investigate if steric 
interaction arising from molecular shape alone can 
induce tilt in the smectic phase. 

The model molecule is comprised of three rigidly 
linked hard sphcrocy finders of length/breadth ra- 
tio L ror x/D arranged in a zigzag configuration (see 
the inset in figure 1). Both ends of the molecule is 
making an angle \& with the core. The idea behind 
using such hard core model is that liquid crystal 
phase behavior is largely entropy driven and de- 
termined by the hard core repulsion between the 
liquid crystal mesogens. Hard spherocylinder pro- 
vides a simple model both in terms of computa- 
tional ease and theoretical approach. It has been 
studied extensively and exhibit a rich phase be- 
havior including isotropic, ncmatic, smectic A, and 
crystal phases. The advantage of using such hard 
core model is that we can vary the shape of the 
molecule (by changing ^) and see how that affects 
the large-scale organization of the liquid crystal 
phases. 

For convenience we introduce reduced units. 
The reduced pressure P* is defined as P* = 
(3Pv] lsc and a reduced density p* = pvhsc, where 
Vhsc is the volume of the straight hard spherocylin- 
der of length L and breadth D. We have per- 
formed MC simulation in NPT ensemble with pe- 
riodic boundary condition on a system of 400 zig- 
zag shaped molecules. The simulation cell consists 
of N = 400 molecules in a cubic box of dimen- 
sion L x x L y x L z . Initially the system was pre- 
pared in a crystal phase at high pressure. Among 
the possible different crystalline order, we choose 
the antipolar crystal ordering corresponding to the 
highest packing density to ensure the highest sta- 
bility. Starting from crystal phase at high pressure 
(P* = 13) we decrease the pressure successively 
by steps of AP* = 1, until we reach a reduced 
pressure P* = 1. For each run, at a given pres- 



from the previous higher pressure is used as the 
initial state. At each state point (P*, V?) the sys- 
tem is equilibrated for 200000 MC cycles and one 
million MC cycles are used for the production of 
the results and the analysis of the various thermo- 
dynamical and structural quantities. During each 
MC step each molecules were chosen randomly and 
displaced using Metropolis criteria. The reorienta- 
tion move was performed using quternion j2^. In 
all the simulations reported below, we use a length 
to breadth ratio L\> an /D = 2. 

In order to fully characterize different phases of 
the systems various order parameters were com- 
puted. The location of the solid-liquid phase 
boundary is determined by computing the squarcd- 
magnitude of the in-layer translation order param- 
eters pGk defined as 



PGk = 



1 M 

— J^expCiGfcTj) 
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FIG. 1: Phase diagram of zigzag shape molecules with aspect 
ratio L/D = 2 as a function of opening angle $ and reduced 
pressure P*. The following phase are present: isotropic liq- 
uid (I), nematic (N), smectic A (SmA), smectic C (SmC), 
columnar (Col), tilted crystal (XT), crystal (X). 



Gi, G2, G3 are the reciprocal basis vectors and 
Yj is the position of the center of mass of the 
molecule j and M is the number of molecules in a 
given layer. 

The smectic-nematic phase boundary is deter- 
mined by the squared-magnitude of the layer trans- 
lational order parameter p^ defined as 



P\\ 



1 N 

— 5^exp(iG|| Tj 



(2) 



j'=i 



where Gn = =j- z, the layer normal being along 

z. 

To distinguish between a tilted and a non tilted 
phase (smectic or crystal) we introduce the in-layer 
polar order parameter m defined as 



m = 



1 M 

— Y 
m f-r 



m ; 



(3) 



where rhj is the unit vector contained in the plane 
of the molecule and passing through one of the 
apex of the molecule. With this definition, care 
must be taken to ensure that two zigzag molecules 
having the z-component of the end-to-end vector 
pointing in opposite direction have in fact the same 
polar direction. 

The orientational order-isotropic phase bound- 
ary is determined by the eigenvalues of the second- 
rank tensorial orientational order parameter Q a p 
defined as 



N ^ \ 2 



(4) 



with a, (3 = x, y, z and n., is the molecular end-to- 
end unit vector of molecule j. The nematic order 
parameter S is given by the largest eigenvalue of 



zero in the isotropic phase and will tend to one in 
highly ordered phase. 

To distinguish between SmA and SmC we com- 
pute the tilt angle of the central spherocylinder 
(core) of the zigzag shaped molecule with the layer 
normal. 

To map out the complete phase diagram we have 
performed simulations for several zig-zag angles 
in between \f = and 80°. The phase diagram 
(\I>, P*) is presented in Fig. [I] 

The system exhibits rich phase behavior with 
isotropic liquid (I), nematic (N), smectic A (SmA), 
smectic C (SmC), columnar (Col), tilted crystal 
(XT), crystal (X) phases. Configurations from 
the isotropic, nematic, smectic A and tilted crys- 
tal phases are shown in Fig. [21 for a zigzag an- 
gle VP = 15°. Also shown on the same figure 
a configuration from the smectic C phase for a 
zigzag angle = 65°. The phase diagram is com- 
pletely isomorphous to the phase diagram obtained 
for bent-core molecules represented as hard-core 
dimer formed by two interdigitated hard-core sphe- 
rocylinders sharing one spherical end cap f 2 ll | . The 
SmC replaces the polar SmA (SmAP) phase found 
in bent-core molecules since polar symmetry break- 
ing leads, for the zigzag model, to a phase having 
all the symmetries defining a tilted smectic phase. 
Close-packing interactions combined to the molec- 
ular zigzag geometrical shape induce a tilt of the 
molecules with respect to the layer normal. 

The nematic phase is stable for zigzag angles 
smaller than <3> = 40°. With increasing zigzag an- 
gle, the region of stability of the nematic phase 
decreases, vanishing for zigzag angles larger than 
~ 40°, leading to an (I, N, SmC) triple point near 
$ = 40°. The existence of a biaxial nematic phase 
remains an elusive possibility in thermotropic LCs. 
Zigzag molecules are good candidate to exhibit 
biaxial nematic due to their geometrical phase. 
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FIG. 2: Final configurations from Monte Carlo simulations of 
N = 400 zigzag molecules as a function of pressure. From left 
to right, top to bottom: crystal (ty = 15°, P* = 11), smectic 
A (* = 15°, P* = 9), smectic C (* = 65°, P* = 8), 
nematic phase = 15°, P* = 5), isotropic phase = 
15°, P* = 1). 

study, confirming the results obtained with bent- 
core molecules . As is clear from simulations of 
the hard biaxial ellipsoid system [2f| [2(| , a biaxial 
nematic phase requires a highly biaxial molecular 
shape. Because such a requirement could be satis- 
fied for larger L/D and specific zigzag angle in the 
present model, the possibility of the presence of a 
biaxial phase should not be ruled out. 

Because straight sphcrocylinders do not exhibit 
any SmC ordering, it is expected that our model 
should exhibit a transition from SmA to SmC. This 
transition occurs for a zigzag angle between 26° 
and 28°, and is associated with two triple points, 
a (SmA, SmC, N) triple point near ^ = 28° and 
a (SmA, SmC, XT) triple point near f = 26°. 
The SmA-SmC transition is purely entropy driven 
and arise from packing effect. In figure 01 we have 
plotted the smectic tilt angle (defined as the av- 
erage over the tilt angle with respect to the layer 
normal made by the central spherocylinder of ev- 
ery molecules) as a function of the zigzag angle. 
This is the first time that an idealized model sys- 
tem gives conclusive evidence of the occurrence of 
a SmC phase. Due to the weak coupling between 
adjacent layers in the SmC phase, it was impossi- 
ble to determine the relative stability of synclinic 



FIG. 3: Evolution of the squared-magnitude of order param- 
eters as a function of reduced pressure for an opening angles 
of *P — 20° (top) and $ = 65° (bottom) showing, respec- 
tively the phase sequences XT-SmA-N-I and XT-SmC-I as 
a function of decreasing pressure. The following order pa- 
rameters are plotted: (A, V, <) solid-liquid order parameters 
Pi, pi, pa', (□) polar order parameter m; (O) smectic order 
parameter pu ; (O) the largest eigenvalue of the nematic order 
parameter Q a /3- 

order by direct simulations. But it is reasonable to 
assume that the zigzag model obeys the same gen- 
eral thermodynamic mechanism like the hard sphe- 
rocylinder and the hard-dimer models that we have 
recently studied 0, 01 ■ I n these studies within a 
framework of "sawtooth" model, we have demon- 
strated that the entropy content of molecular-scale 
fluctuations of the interface between smectic C 
layers ("out-of-layer" molecular fluctuations) pro- 
vides a general thermodynamic mechanism that 
uniquely favors synclinic ordering. 

We also find a transition between a tilted crys- 
talline phase and a untilted crystalline phase (i.c 
a rotator phase). This rotator phase is stable for 
zigzag angles smaller than ^ = 20°, and is char- 
acterized by a (SmA, X, XT) triple point around 

= 18°. Quite interestingly, the rotator phase 
competes with a columnar phase for zigzag angles 
larger than 'J ~ 3° and smaller than \I> ~6°. This 
narrow columnar phase is characterized by signifi- 
cant two-dimensional crystal order parameters but 
a negligible magnitude of the smectic order param- 
eter. 

Insights into the shape of the phase boundaries 
can be gained by supposing, to a first approxima- 
tion, that the partition function of the system can 
be decomposed into a product of positional and 
orientational contributions, in which case the en- 
tropy is the sum of an orientational entropy and 
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FIG. 4: Tilt angle (0) (top) and layer spacing (bottom) as a 
function of zigzag angle. 



ferent forms of entropy determines the stability of 
a given phase at a given density. In the limit of 
straight spherocylinders = 0°), the isotropic- 
nematic phase transition occurs when the gain in 
positional entropy S pos exceeds the loss of orien- 
tational entropy S" 31 ' 1011 * (27}. A nematic-smectic 
phase transition occurs when the gain in transla- 
tional entropy perpendicular to the long molecu- 
lar axis exceeds the loss of positional entropy 
parallel to the long molecular axis S pos , leading 
to the formation of a stack of two-dimensional liq- 
uid layers. Similar reasoning can be applied to 
zig-zag molecules: in the range 0° < \& < 40°, 
the isotropic phase is more favourable at smaller 
zigzag angles. As the cores become more bent 
(larger zigzag angles), the gain in positional en- 
tropy associated with nematic ordering is reduced, 
and the nematic phase range is reduced, eventu- 
ally disappearing for \I/ > 40°. The shape of the 
nematic-SmC boundary (i.e., for 30° < * < 40°) 
can be qualitatively understood in the same way 
by noticing that the positional entropy parallel to 
the long molecular axis S^ OB is larger for smaller 
zigzag angles than for larger ones, stabilizing the 
nematic phase for smaller zigzag angles. 

The model exhibits a rich phase behavior includ- 
ing tilted and nontilted crystal, columnar, smectic 
A, smectic C, nematic and isotropic phases. The 
model shows without any ambiguity that excluded 
volume interaction arising out of molecular shape 
is sufficient to produce tilted smectic phase even in 
the absence of electrostatic interactions. 
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